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Abstract. - We discuss the optical loading of a Bose-Einstein condensate in the Thomas- Fermi 
regime. The condensate is loaded via spontaneous emission from a reservoir of excited- state 
atoms. By means of a master equation formalism, we discuss the modification of the condensate 
temperature during the loading. We identify the threshold temperature, Tth, above (below) 
which the loading process leads to cooling (heating), respectively. The consequences of our 
analysis for the continuous loading of an atom laser are discussed. 



The experimental realization of a Bose-Einstein condensate (BEC) has aroused an enor- 
mous interest during the recent years pp. Among the disparate variety of topics related with 
the BEC physics, the possibility of accomplishing a coherent source of atoms, or atom laser, 
has attracted a large attention As a coherent source of matter waves, the atom laser is 
expected to lead to new applications in atom optics, and its relevance is often compared to 
the impact of optical lasers in standard optics. The atom-laser has been produced experimen- 
tally after outcoupling atoms from a BEC either by employing radio-frequency fields 00], 
or Raman techniques [5] . The coherent character of the source has been demonstrated in a 
series of interferometric experiments 0JEj][7|. The availability of continuous- wave (cw) atom 
lasers would open the way to "high power" and precision applications. Several groups have 
developed various techniques to outcouple atoms in a (quasi-) continuous way However, 
the continuous outcoupling represents just a half way towards a cw atom laser. Without a 
continuous refilling of the BEC, the atom laser output lasts only as long as some atoms in 
the BEC are kept. In addition, the continuous refilling of the condensate could be employed 
to repair a condensate against inelastic losses, and in this way it could allow to analyze the 
BEC physics at larger densities, at much longer time scales. Recently, the first steps toward 
this direction has been achieved in an experiment performed at MIT [H]. 

Two different physical mechanisms could provide a continuous pumping into a BEC. On 
one hand, collisional mechanisms [5], in which two non-condensed reservoir atoms collide, one 
being pumped into the BEC, whereas the other carries most of the energy and is evaporated. 
Several proposals have been suggested in this direction. Ref. analyzes the continuous 
evaporative cooling of an atomic beam directed along a long magnetic guide, at the end of 
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which the atoms are predicted to enter the degenerate regime. This method demands a large 
flux of slow atoms in order to achieve a sufficiently large collisional rate which guarantees the 
required efficiency of the evaporative cooling in a finite guide. Recently, Roos et al. jlOj have 
analyzed a scheme in which particles from an incident beam are trapped in a potential well 
when colliding with particles already present in the well. For a cigar-shape potential, it has 
been shown that a resonance can occur when the transverse evaporation threshold coincides 
with the energy of the incident particles, leading to a large increase in phase space density with 
respect to the incident beam. We should also refer at this point to a recent experiment f||, 
where, by means of an optical tweezer, a BEC is merged with a pre-existing one. A second 
way of continuous loading of a BEC is provided by the spontaneous emission from a reservoir 
of excited-state atoms |llU12lll3] . If this reservoir could be filled in a (quasi-) continuous way 
by laser cooling techniques, one would benefit from the large cooling efficiency of laser cooling 
compared to evaporative cooling, allowing for a considerable increase in atomic flux produced 
by an atom laser. In order to achieve the continuous optical loading of a BEC, it is crucial 
to avoid the heating introduced by the photon reabsorption [14) . Fortunately, it has been 
shown that in the so-called Festina Lente regime, in which the spontaneous emission rate, 7, 
is smaller than the frequency of the harmonic trapping, u>, the heating due to reabsorption is 
practically suppressed jl5j . It has been also recently reported that, under certain conditions, 
the reabsorption can be also largely suppressed in atoms with an accessible three level A 
scheme ^21- Additionally, quantum interference effects can lead to positive (cooling) effects 
of the reabsorption in the so-called BAR regime |16II17| . 

The continuous optical loading of a BEC have been previously analyzed in ref. |12| , in which 
the BEC was considered in the weak condensation regime, where the mean interaction energy 
is smaller than the separation between the harmonic levels. In that case, the loading dynamics 
can be split from the collisional physics in the corresponding quantum Master equation (QME), 
and both phenomena can be simultaneously simulated via Monte Carlo methods 18, 12 . In 
particular in ref. |12j , it was shown that a BEC can be first created in an initially empty trap 
by means of optical loading, and then maintained against losses by properly fixing the pumping 
rate. However, the weakly-interacting regime sets important restrictions to the atomic density. 
In the present paper, we extend our analysis to arbitrary regimes, and in particular to the 
Thomas- Fermi regime, in which the mean interaction energy is much larger than the separation 
between the trap levels. In that case, the BEC wavefunction and the excitations depend on the 
particle number, and as consequence the wavefunctions vary dynamically during the pumping 
process, complicating the analysis of the problem. In this paper we develop the necessary 
theory to analyze the problem. In particular, we study the modification of the condensate 
temperature during the pumping process. We show that there is a threshold temperature for 
a given number of trapped particles above (below) which the system is cooled (heated), and 
show that this effect could be employed to control the condensate temperature. 

We consider a sample of cold atoms of mass m with an accessible electronic two-level 
system, formed by the ground state \g) and an excited state |e). Both \g) and |e) atoms are 
in a harmonic trap, which for simplicity is considered isotropic, with frequencies u> g and uj e , 
respectively. The \g) atoms are Bose condensed, i.e. with a temperature T <C T c , where T c is 
the critical temperature for onset of the condensation. The formation of the condensate using 
spontaneous emission was previously discussed in ref. ^21- The |e) atoms, which eventually 
decay spontaneously into the \g) state, are considered as thermally distributed. 

The atoms in the \g) trap are described by the corresponding field operator $ s (r). For 
sufficiently low T, when the number of condensed atoms in the \g) trap, No, is comparable to 
the total number of atoms, A, one can employ the Bogoliubov approximation and substitute 
i& g by a c- number, ipo, the BEC wave function, whose stationary value is provided by the 
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time-independent Gross-Pitaevskii equation (GPE) 

HgpM?) = (-|^ v2 + + U gg N\^ \ 2 - m) = 0, (1) 

where p is the chemical potential, and U gg = 4TTh 2 a gg /m is the coupling constant for the 
collisions between \g) atoms, with a gg the scattering length. In eq. we have neglected the 
e-g collisions, since the number of atoms in the |e) state is assumed very small. 

The excitation spectrum is obtained after linearizing i& g around the ground-state solution, 
5' g (r) = ipo(r) + Sjp(f f ). The perturbation Sip is then expanded in the standard form |19| 
5ip(r) = Ufiffigfi — Vft(r)gt, where g^ and gt are the annihilation and creation operators for 
the quasiparticles with spherical quantum numbers n — n,l,m. The wave functions Uff(f) 
and Vfi{r) obey the standard orthonormalization rules [19] ■ Linearizing in the ua(r) and Ujj(r) 
wavefunctions, one obtains the corresponding Bogoliubov equations 

{H GP + U gg N Q Qi, 2 Q)u n + U gg N Qip 2 Q*v n = fiCjiun (2) 
-U gg N Q*Mr) 2 Qun - (H GP + U gg N QiplQ)*Vfi = hCj 9 n v H (3) 

where we have used the projection operators Q = 1 — \ipo)(ipo\ (Q*) orthogonally to ipo 
(V'o) E2> an( i &n denotes the quasiparticle excitation spectrum. 

The physics of the atoms in the |e) state is described by the Schrodinger equation: 

n^(f) = (~|^V 2 + |c e V + 2U ge N \^ tfn(r) = m^kir), (4) 

where we have taken into account that the |e) atoms are affected by the mean- field potential 
induced by the collisions with the \g) atoms, which are characterized by a coupling constant 
U ge = Anh 2 a ge /m, with a ge the corresponding scattering length. Due to the low density in 
the |e) trap, we have neglected in eq. (@J the e-e collisions. In the following, et and denote 
the creation and annihilation operators in the eigenstate ip^, with eigenfrequency u)^. 
The interaction of the atoms with the vacuum electromagnetic field is given by 

rh V \ n n / 

where at (a^) are creation (annihilation) operators of photons of wavenumber k, polarization 
is, and frequency u> u {k). In eq. (J5J, d is the atomic dipole, t^ v is the polarization vector, and 
p(k) the density of states. The Frank-Condon (FC) factors = J dfip^(r) exp(ikr)ip (f) 
characterize the transitions into the condensate, whereas 77-^ — J dfip^r) exp(ikr)Uft(r) de- 
termine the transition into the u quasiparticle wavefunctions. In the same way, the coefficients 
r)%-t characterize the transition into the v wavefunctions. Finally, the vacuum energy is pro- 
vided by the Hamiltonian: TCf = J~) v J dkuj u (k) at a% u . The physics of the system is therefore 
described by the Hamiltonian TL — TL a + H a f + rlf, where the atomic Hamiltonian can be 
written in the form (extracting the condensate energy) H a = Yla ^5^5n+Sm( a; o+^)et^e^ 
where ojq is the transition frequency between \g) and |e). 

From the previous Hamiltonian, and using standard techniques of quantum stochastic 
processes [20] ^ we get in Born-Markov approximation the QME for the density matrix p: 



p(t) = -iH eff p + ipHljf + Jp, 



(6) 
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where H e ff is the effective (non-Hermitian) Hamiltonian, and 

I m m,n rn,n J 

is the jump super operator, where ^ niV rh'(u ) = I dnw {^)v%,n{^k)v^ n' ( w oe fe )* (a=c, u, 
v), with W(f2) denoting the emission pattern, and 7 is the spontaneous emission rate. 

The effective Hamiltonian H e ff is related with photon reabsorption processes, whereas the 
jump superoperator describes a spontaneous emission act without any further reabsorption. 
In this paper, we consider the reabsorption processes as negligible. As commented above, this 
should be the case if the Festina Lente condition is fulfilled [2 , or alternatively for atoms with 
a three-level A scheme under the conditions of ref. (T2|. Therefore, we constrain ourselves to 
the analysis of the jump processes from |e) to \g). From eq. (J7J) the amplitude for a jump from 
a state n e of the |e) trap to a state fig (a = u, w, c) of the \g) one is given by 

r 3 (.j He = 7/ dnW(mn { g a) \e^\rie)\ 2 , (8) 

with (r\n e ) = R le n e (r)Y lm (n), and (r\n g a) ) = (r)Y Jm (fi). The FC factors are of the form 

OO i 

l"=0 m=-l" 

where -^J, j"„ e j e = / rf^ 2 Pi g n g ( r ) 3i"(kr)Ri cnc (r) is the radial integral. The angular part is 
provided by Bj" m = J c2Q r Y5* m (O r )Y/" TO (Q r )F; em<! (Q r ), solvable in terms of the corre- 
sponding 3-j symbols. Due to the properties of these symbols, Bj m lenie is independent of m 
since — m g + m e + m = 0, l g + l e + I" = 2p (where p is an integer), and \l g — l' e \ < I" < l g + l e - 
From the amplitudes JHJ and the populations of the corresponding levels, we obtain the tran- 
sition probabilities describing creation of a quasiparticle n — (n g ,l g ) 

PZ = ((N n (T)) + l)Y,(Nfi e (Te)) £ r 4"»,n E - ( 10 ) 
n e i c m g ,m c 

Similar expressions are obtained for the transition probability of destruction of a quasiparticle 
n associated with the transfer of a particle into the condensate after thermalization, P%, and 
the direct decay process into the condensate, P c . In eq. (fTUfl . (N n (T)) = (exp[E n /K B T] — 
is the Bose-Einstein distribution in the \g) trap, and (iVg e (T e )) = N exp[—Efi e / KbT b ] is the 
Boltzmann distribution (normalized to 1) in the |e) trap. 

Knowing the transition probabilities P% (|10f) . P% and P c we can simulate the effects of 
a spontaneous emission process. We assume that the time scale of thermalization in the \g) 
trap is much faster than the time interval between different pump events, i.e. the collisional 
rate in the ground-state trap is much larger than 7. In this way, after every pump, the system 
re-thermalizes, redistributing the energy gained or lost during the pump process. We analyze 
in the following the variation of the condensate temperature after the pumping. 

After pumping, the number of \g) atoms increases to N' = N + 1, and the energy becomes: 



(E') =E + M P C + 5>„ + (iC n (P: - o 

n 



(11) 
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where E is the energy before the pumping, and C n = |ti n | +|t> n | . After the re-thermalization, 
the system acquires a new temperature T' which can be evaluated from the expressions: 

E' = (^0|H G -^|^0| 2 |V'0}+E Ar «( T ')( e »+^n)-E( e «-/ i )l W ™| 2 ' ( 12 ) 

n n 

N (T') = A + 1-]T7V„(T')C„ + M 2 , (13) 

n 

where N n {T) = (exp[e n /KBT] — 1) _1 is the Bose-Einstein distribution in the \g) trap. In 1)12(1 
and (|13|l . we assume the excitation spectrum as that calculated for N particles, i.e. before the 
pumping. This approximation is valid as long as the number of excited particles N — Nq ;§> 1. 

Assuming T' = T + ST with ST <C T, we can perform a Taylor expansion up to first 
order in ST/T. Substituting into eqs. 1(1 2() and 1(13(1 . we can obtain the expression for £ = 
(Nq/N' — No/N)/(Nq/N). If £ > 0, the pumping increases the relative population in the 
condensate, and therefore the temperature of the sample is reduced. On the contrary, if £ < 
the system is heated by the pumping. It is easy to see that up to terms of the order SN/Nq, 
where 8N — N — No, the system is cooled (£ > 0) when the condition: 

( N ^ &W- pc )-Y, n ^ + jiC n ){PZ-PZ) sr exp[e n ] 

( ° " } y V- c *p^>' , ^ £ " C "(exp[ e ~„]-l)2 > 5N (14) 

is fulfilled, with jl — (i/ KbT , e n = e n / KbT . For moderately low T this condition is typically 
fulfilled. For large T, SN can increase enough to violate 1(141) . On the other hand for very low 
T the transition probability into the condensate P c — ► 1, and also in this case the conden- 
sate is heated. Hence, there is a transition from heating to cooling at some finite threshold 
temperature, Tth, which depends on four different parameters: the total number of atoms N, 
the temperature, T e , of the excited state trap, the coupling constant U eg for the collisions 
between |e) and \g) atoms, and the ratio oj e /ujg between the frequencies of the \g) and |e) 
traps. When N increases, so does jl and it is easy to see that T t h decreases. If T e increases 
the transition to higher lying excited states are more favorable, and as one could expect T t h 
increases. The parameters U eg and uj e /ujg influence the spectrum of the excited state trap. 
For U eg — the harmonic oscillator spectrum is recovered. However, when U eg grows, the 
mean- field repulsion leads to an expulsion of the |e) atoms from the trap center, increasing the 
probability for a transition into excited states of the \g) trap. On the other hand, if uj e 3> u) g , 
and T e is sufficiently small, the mean-field shift just provides an overall homogeneous shift of 
the levels of the |e) trap, and the T t h becomes comparable to that for U eg = 0. 

We have analyzed ST after a pumping using the following procedure. For given values of 
TV and T, we first calculate the number of atoms in the BEC, Nq, and evaluate ipo by evolving 
eq. Q in imaginary time. Then, we diagonalize eqs. © and JSJl employing a harmonic 
oscillator basis, to find the {u, v} eigenf unctions, and the corresponding eigenenergies ( : ) . 
As a next step, we calculate r_( ) - (a — c,u,v) using eq. iJSJ. For a given T e we evaluate 

P^, and P c . From eq. ((11(1 we obtain the new average energy after the pumping. Finally, 
from eqs. 1(12(1 and 1(13(1 we obtain the new temperature after the rethermalization. 

Fig. nji shows the variation of T/T c after a spontaneous emission as a function of the 
initial T/T c for lo = 20 x 2tt Hz, u e = 10cj s , T e = 2lfiuj gi and N = 10 4 , 5 x 10 4 and 
10 5 . In the calculation of the FC factors 10, we have considered a Lamb-Dicke parameter 

( 1 )In TF regime u n ,v n and e n are known analytically. We employ here however the numerical approach, 
which works in all intermediate regime between weakly condensed and TF regime 
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NxlO 4 NxlO'' 



Fig. 1 - a) variation of T/T c after a spontaneous emission as a function of the initial T/T c for 
uj = 20 x 2tt Hz, uj e = lQu g , T e = 2lhcv g , and N = 10 4 , 5 x 10 4 and 10 5 ; b) dependence of T th on 
N for u e — uj g and U eg /U gg = l(dotted-dashed curve), 0.5 (dotted curve), and (dashed curve); 
the solid curve shows the case u> e = WuJ g and U eg /U gg = 1; c) solid line indicates the Monte Carlo 
simulation of the loading process and dashed line T t h for uj e = 10oj 9 , T e = 21huj g 



i] 2 = E rec /hu>g = 1, where E rec is the photon recoil energy. The choice of 77 and T e is due 
to our numerical limitations, but we expect qualitatively similar results for other values of 
these parameters. As previously discussed, the system is cooled for T > T t h, where T t h is 
progressively lower for larger N. In fig.^i, due to the low value of T e chosen, it is not possible 
to observe the heating region for large temperatures, which for this case is located at T close 
to T c . For other ranges of parameters the qualitative picture does not change significantly, 
i.e. always a threshold temperature separating heating and cooling regions exists. 

We have analyzed the dependence of T t h on different parameters. Fig. ^ shows the 
dependence of T t h on N for to e — to g and different values of U eg . For U eg — > 0, T t u decays 
monotonically with TV. However, when U eg increases, and due to the repulsive mean-field 
induced by the e-g collisions, the eigenfunctions of the |e) trap possess a local minimum at 
the trap center, and the pumping into the condensate becomes less effective. Consequently, 
T t h increases for large N, as observed for U eg — U gg . This effect is less pronounced, for larger 
values of u e keeping fixed T e , since then the lower levels of the |e) trap have an extension 
smaller than the BEC wavefunction, and hence the mean-field just produces a global energy 
shift of the lowest |e) levels, with no consequences for the pumping (see fig. OJj). 

Finally, we have simulated for different conditions the continuous optical pumping into the 
BEC. To this aim, we have employed Monte Carlo methods to evaluate the corresponding rate 
equations, with transition probabilities (f 101) . P T " ; and P c . In principle, after every pumping 
the previously described algorithm should be repeated. In practice, it is enough to do so after 
a number of pumps N pumps -C 6N, which is dynamically adjusted in our calculations. Fig. 
shows our results for u e = 10w s , T e — 1\hio g . The dashed curve in fig.^ indicates T t h- The 
simulation is started within the heating region (below the T t h curve), with TV = 4.8 x 10 4 , and 
T/T c — 0.18. As expected, T increases with the pumping. However, once (N,T/T C ) crosses 
the T t h curve, the temperature begins to decrease when pumping. Therefore, interestingly, 
the temperature of the system should asymptotically approach the curve T t h- 

In this paper, we have analyzed the continuous optical loading of a BEC in the Thomas- 
Fermi regime. Contrary to the weak-condensation case, the BEC and excited-state wave- 
functions vary during the loading. By means of GPE and Bogoliubov equations, we have 
determined before every pumping the proper wavefunctions, and transition probabilities. As- 
suming rapid thermalization, we have monitored the variation of the temperature during the 
loading. We have observed that for a given number of trapped atoms, N, there is always a 
threshold temperature, T t h, below which the BEC is heated. The concept of T t h is very general 
and should also apply to collisional loading of a BEC. We have analyzed the dependence of T t h 
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on different parameters, in particular N, the trap geometries, and the interparticlc interac- 
tions. A lower T t n is obtained for smaller interactions between excited-state and ground-state 
atoms, since repulsive interactions displace the excited-state atoms away from the trap center, 
preventing the effective pumping into the BEC. Our analysis shows that, not only the number 
of atoms in the trap, but also the temperature of the system, can be maintained by means 
of optical loading. To this aim, the condensate should be created with a temperature and a 
number of atoms, lying close to the curve T t h{N). In that case, in the presence of losses, if the 
loss mechanism leads to heating, T becomes larger than T t h(N), and the system is cooled by 
next pumpings, whereas the opposite is true if the loss mechanism cools. Interestingly, this 
could allow for the sustained analysis of BEC at a quasi- const ant temperature. 

Several physical scenarios could be devised for the continuous optical loading of a BEC. 
For example, an optical lattice could be used to move atoms in an internal state |r) from 
a relatively hot reservoir into the center of the BEC region, where a Raman pulse could be 
employed to transfer the \r) atoms into an state |e), which decays into the ground-state \g). In 
this way, it could be achieved that only those lattice sites near the condensate center actually 
decay into the \g) trap. For a sufficiently strong optical lattice each lattice site will behave as 
a single trap with few occupied levels, as that discussed in our paper. Interestingly, in this 
scenario, due to the small size of the |e) trap, T t h should monotonically decay with N, since 
the mean- field effects should just shift globally the states of the \e) trap. Therefore, in this 
system the pumping could be employed to post-cool an already formed BEC down to T = 0, 
and maintain the condensate temperature very low against possible loss sources. 

We acknowledge support from the Alexander von Humboldt Stiftung, the Deutsche Forschungs- 
gemeinschaft, EU RTN Network "Cold Quantum Gases", ESF PESC Program BEC2000+. 
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